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^ I Abstract 

■ We calculated profiles of CIV 1550, Si IV 1400, NV1240 and OVI 1035 doublet 






lines using results of 3D MHD simulations of disc accretion onto young stars with 
dipole magnetic field. It appeared that our calculations can not reproduce profiles 
of these lines observed (HST/GHRS-STIS and FUSE) in CTTSs's spectra. We also 
found that the theory predicts much larger C IV 1550 line flux than observed (up to 
^^. two orders of magnitude in some cases) and argue that the main portion of accretion 



energy in CTTSs is liberated outside accretion shock. We conclude that the reason 
H ' of disagreement between the theory and observation is strongly non-dipole character 

^ . of CTTS's magnetic field near its surface, a 



o 
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1 Introduction 



> 

^ ■ Since the beginning of the 1990s, there has been a consensus that the hne and continuum 

^ . emission observed in the spectra of classical T Tauri stars (CTTSs) results from the 

magnetospheric accretion of circumstellar material. More precisely, the magnetic field of 
the star is believed to stop the accretion disk from reaching the stellar surface. In some 
way the disk material becomes frozen in the magnetospheric field lines and slides along 
them toward the stellar surface, eventually being accelerated to velocities Vq ~ 300 km/s. 
The gas is then decelerated in an accretion shock (AS), whose radiation presumbly gives 
rise to the observed line and continuum emission. 
5^ I Radial extention of pre- and post-shock radiating regions of CTTS's AS is much smaller 

than stellar radius - Zpre, Zpst -C -R* in figure [H - making it possible to calculate the 
structure and spectrum of the AS in 1-D approximation (Lamzin, 1995). Calculations of 
Lamzin (1998) and Calvet & GuUbring (1998) indicated that the structure of the flow can 
be specifled nearly unambiguously by two parameters: the velocity Vq and density po (or 
particle number density A^o) of the gas far in front of the shock. 



^ Published in Proc. of lAU Symp. 243: "Star-Disk Interaction in Young Stars", Grenoble, France, 
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Figure 1: Schematic structure of CTTS's AS. 



Calvet & Gullbring (1998) used results of their calculations to derive the parameters 
of the accretion shock via modeling of the continuum spectral energy distributions of 
classical T Tauri stars. However, they did not took into account limb darkening effect 
and what is more the agreement between the calculated and observed spectra of the veiling 
continuum cannot be considered as a decisive support for the magnetospheric model, since 
boundary-layer models provide equally good agreement - see e.g. Basri & Bertout (1989). 
The line spectrum is far more informative, and comparisons of the calculated and observed 
intensities and profiles of emission lines enable detailed studies of the accretion processes. 

Optically thin lines are best suited for this purpose: their intensity ratios can be used 
to derive physical conditions independent of the geometry of the region where they are 
formed, and the line profiles provide information about both the velocity field and geom- 
etry of this region. The calculations of Lamzin & Gomez de Castro (1999) demonstrated 
that the O III] 1663, Si III] 1892 and C III] 1909 lines should display the highest intensities 
among the optically thin lines, and just these lines were used to determine the accretion- 
shock parameters for several young stars in that paper. However, Gomez de Castro & 
Verdugo (2001, 2007) questioned whether these lines form in the AS. Spectral lines of 
neutral or singly ionized atoms apparantly form not only in AS but in magnetospheric 
flow and wind as well - see e.g. Edwards (2007) - and thus likewise cannot be used for 
diagnostic of AS. 

Resonant UV lines of the C IV, Si IV, N V and O VI uvl doublets looks as the most 
suitable lines for diagnostic of CTTSs's AS, especially lines of CIV 1550 doublet: they 
are strong in CTTS's spectra and calculation of their intensities is relatively simple. 




Figure 2: Relative intensities of CIV 1550 doublet lines expressed in %. See text for 
details. 



2 Intensities of C IV 1550 doublet lines: theory vs. 
observations 

Generally speaking, lines of C"*"^, Si"*"^, N"*"^ and O"*"^ ions form before and behind the shock 
front. Gas temperature in the pre-shock (precursor) zone of CTTSs does not exceed 20.000 
K, but ions up to O"*"^ (at 200 < Vq < 400 km/s) exist here due to photoionization of 
accreted matter by X-ray and UV quanta from post-shock cooling zone. When infalling 
gas crosses the shock front its temperature raises up to 1-3 MK and for example C"*"^ ions 
almost complitely transform to C^^ ions. Then gas cools and ions of interest appeares 
again but at that moment gas velocity is very close to zero - see Lamzin (1998) for details. 

Lamzin (2003a) carried out non-LTE calculations of profiles of CIV 1550, Si IV 1400, 
NV 1240 and O VI 1035 doublet lines for a plane-parallel shock viewed at various angles. 
Calculations were performed for the range of preshock gas parameters Vq, Nq appropriate 
for CTTSs. Intensities of CIV 1548+1551 lines, normalized to JF = pqVq/A value, as a 
function of cosine /i of an angle between the normal to the shock front and the line of 
site are presented in figure El The ratio is expressed in percent, such as different lines in 
fugure respect to different infall gas velocities (Vq = 200, 300 and 400 km/s). Presented 
results were calculated for gas particle density A^o = 10^^ cm~^ (left panel) and 10^^ cm^^ 
(right panel). 

The value JF was chosen for normalisation because just this value expected to be equal 
to the bolometric flux of veiling continuum emission produced by AS. Indeed, according 
to current paradigm a half of X-ray and UV quanta from post-shock cooling zone moves 
to the stellar surface. These quanta are absorbed in upper layers of stellar atmosphere 
and then should be reradiated predominantly in the continuum - see Calvet & Gullbring 
(1998) for details. Thus the ratio S, depicted in figure [21 is theoretical prediction for ratio 
of C IV 1550 doublet line flux to bolometric flux of veiling continuum. Thus this ratio 



expected to be ~ 1 % almost independenly on Vq, A^o parameters of AS. 

Meanwile it was found that observed ratio is much smaller - see table [T] in which we 
summarised results of our analyses of CTTSs UV spectra observed from Hubble Space 
Telescope (Kravtsova & Lamzin, 2002a,b; Kravtsova, 2003; Lamzin et al., 2004). The 
descripancy between the theory and observation is significant - more than two orders 
of magnitude in the case of RY Tau and DR Tau, such as this conclusion does not 
depend on current uncertainty of value and law of interstellar extinction in the direction 
to investigated stars. 





Table 1: Relative contribution of CIV lines in emission of AS 


Star: 


RY Tau DR Tau T Tau DS Tau BP Tau DG Tau 


Theory 


<5, % 


0.002 0.003 0.02 0.02 0.04 0.07 


~ 1.0 



We suppose that this discrepancy means that the main portion of veiling continuum 
(up to 99 % in some cases juging from table [T]) originates outside C IV 1550 line formation 
region, i.e. outside strong (Mach number M^h ^ 1) AS. In other words we conclude that 
the main portion of accreted matter does not pass throught the AS and falls to the star 
almost parallel to the stellar surface. In this part of accretion flow transformation of 
kinetic energy of infalling gas into the heat and then into radiation should occur in the 
same way(s) as in a boundary layer, i.e in a series of weak {Mgh — 1) oblique shocks. One 
can estimate the maximal possible angle 'jsh between front of the weak oblique shocks and 
stellar surface as follows. 

In the coronal equlibrium approximation C"*"^ ions forms at T ^ 10^ K. Such tempera- 
ture can be reached immediately behind the shock front if accreted gas velocity component 
Vr normal to stellar surface is ^ 70 km/s. Therefore oblique ASs with Vr less than this 
value can not contribute to C IV 1550 line emission of accretion flow but produce veiling 
continuum and emission in lines of neutral, singly or twice ionized atoms. If typical infall 
gas velocity Vq is ~ 300 km/s, then ■jsh = sin~^(K/Vo) < 15°, i.e. gas producing such 
weak shocks indeed falls to the star almost parallel to its surface. 

3 Profiles of C IV 1550 doublet lines: theory vs. ob- 
servations 

As was mentioned above CIV 1550 doublet lines form in two spatially distinct regions of 
strong AS: in the radiative precursor and in the post-shock zone. Gas velocity in these 
regions are different: ^ ~ Vq ~ 300 km/s in the pre-shock zone and ~ 5 — 10 km/s in the 
post-shock line formation region. As a result proflle of e.g. CIV 1548 line in the spectrum 
of plane-parallel shock, viewed from the direction perpendicular to the surface of shock 
front, should have two components: the flrst one is almost at zero-velocity position and 
the second is redshifted to Vq. As follows from our calculations (Lamzin 2003a), both 
components are optically thick, resulting in FWHM of each componet ~ 20 — 30 km/s, 
with relative strength depending on Vq : at Vq < 300 km/s "zero- velocity" component is 
stronger than "high- velocity" one and vice versa at Vq > 300 km/s. If the shock is viewed 



from the direction, that makes an angle 9 with the perpendicular to the shock's surface, 
then "zero- velocity" component should be seen practically at the same position but the 
redshift of "high- velocity" peak should be now Vo/i, where /i = cos^. The same is true 
(in a qualitative way) for lines of Si IV 1400, NV1240 and OVI1035 doublets. 
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Figure 3: Profiles of CIV 1550 doublet lines in spectra of some CTTSs. 



Consider now a part of CTTS's surface occupied with strong AS (accretion zone). Ob- 
served profile of e.g. CIV 1548 line emitted by AS is a sum (an integral) of double-peaked 
profiles from all elementary area AS* of the accretion zone (multiplied to fiAS factor). 
All elementary areas are viewed at different angles due to curvature of stellar surface and 
these angles varies with time due to stellar rotation. One can expect that intensities of 
"zero- velocity" components from all parts of accretion zone will be summarised and the 
(weighted) sum of high-velocity components will results in more or less wide red wing or 
separated redshifted component depending on distribution of Vq, A^o parameters in accre- 
tion zone and on its geometry. Obviously the profile should vary with time due to stellar 
rotation and non- stationary accretion as well. 

Lamzin (2003b) calculated profiles of C IV 1550 doublet lines from strong AS assuming 
that: 1) matter falls to the star in radial direction; 2) Vq and Nq parameters of the shock 
are constant within accretion zone; 3) the zone has the shape of circular spot or sperical 
belt. Results of the calculations were compared with profiles of the lines in UV spectra of 
CTTSs observed with Goddard High Resolution Spectrograph (GHRS) and Space Tele- 
scope Imaging Spectrograph (STIS). Observational data were extracted from Scientific 
Archive of Hubble Space Telescope (http://archive.stsci.edu/hst/target_descriptions.html). 
Calculated profiles differs significantly from observed ones presumbly because our assump- 
tions about the character of accretion flow near stellar surface were not realistic enough. 

One can expect better agreement if to use parameters of accretion flow derived from 
3D MHD simulations of disc accretion to a slowly rotating magnetized young star with 
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Figure 4: Profiles of CIV 1550 and O VI 1035 doublet's components in spectra of T Tau 
(left panel) and TW Hya (right panel). 



its dipole moment inclined at an angle a to the stellar rotation axis. Accretion rate Mac, 
polar magnetic field strength B as well as mass and radius of the central star are free 
parameters of these simulations in addition to the angle < a < 90° - see Romanova et 
al. (2003) for details. Velocity field Vq and gas density po at stellar surface, adopted from 
the simulations, we used as input parameters to calculate profiles of CIV 1550, Si IV 1400, 
NV1240 and OVI1035 doublet lines. For all models we adopted M, = 0.8 Mq, R, = 1.8 
R0, 5 = 1 — 3 kG and varied a, Mac parameters in the 0° — 90° and 10"*^ — 3 ■ 10~^ 
M0/yr intervals respectively. Profiles were calculated for each accretion zone's model with 
different values of an angle between stellar rotation axis and the line of sight (0° < i < 90°) 
as well as for a set of phases of stellar rotation periods, i.e. for different angles ip (in 27r 
units) between magnetic dipole axis and the plane, which contains rotation axis and the 
line of sight (0 <■?/'< 1). 

Observed profiles of CIV 1550 doublet lines in DS Tau, BP Tau, DF Tau and DR 
Tau spectra are shown in figure [3l Solid and dashed lines depicts C IV 1548 and C IV 
1551 componens of the doublet. Profiles of CIV 1550 doublet's components in spectra of 
T Tau are shown in figure H] (left column). T Tau is the only star where there is more 
than one high resolution UV spectrum and one can observe variability of CIV doublet 
lines profiles. Only in the case of TW Hya (right panel of the figure) there is possibility 
to obtain information about lines of O VI 1035 doublet - see Ardila (2007) for reference 
and details. 

We plot in figure[5]results of our calculations for the model with a = 30°, Mac — 4-10^^ 
M0/yr, i = 10° (left panel) and i = 70° (right panel). Vertical row of profiles in each 
panel corresponds to the following set of rotation period phases (from top to bottom): 
ijj = 0, 0.25, 0.5, 0.75. Profiles were normalized to maximal intensity of the line ai ijj = 
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Figure 5: Theoretical profiles of C IV 1548 line calculated for CTTS with dipole magnetic 
field, axis of which inclined at a = 30° angle to the rotation axis of the star. See text for 
details. 



phase. It was assumed that accretion disk does not prevent to observe the part of the star 
that situated below disk's midplane - this is the reason why some profiles have extended 
blue wing. 

Matter falls to the star with dipole magnetic field at the angle 9 < 90° relative to 
its surface. In the absence of magnetic field oblique shock should arise in such situation, 
what means that: 1) the shock front is parallel to stellar surface; 2) velocity component V^ 
which parallel to stellar radius is the pre-shock velocity Vq. But bear in mind that accreted 
gas moves along magnetic field lines it also seems resonable to suppose that shock front 
is perpendicular to the magnetic field lines and therefore Vq = V. To avoid discussion of 
this problem we calculated profiles for both cases: solid lines in figure [5] corresponds to 
profiles calculated for Vq = Vj. case and dashed - for Vq = V case. 

Such approach looks resonable at the moment because both types of theoretical pro- 
files differs from profiles of C IV 1548 line in spectra of CTTS presented in figures El HI 
Observed profiles have only one peak, maximum of which is almost at zero velocity posi- 
tion. The only exception is DR Tau: profile its C IV 1548 line consists of two redshifted 



components but intensity of "high- velocity" component is larger than " low- velocity" one. 
Theoretical profiles calculated for models with another values of a, Mac and i parameters 
have qulitatively the same shape as in figure El i.e. also can not reproduce observations. 
We suppose that the reason of the descripancy is too small divergency of accreted 
gas stream lines within accretion zone that itself occupies only ~ 5 % of stellar surface 
(Romanova et al., 2003). If divergency of the stream (and therefore magnetic field) 
lines within accretion zone would be larger then it seems possible to obtain single-peak 
profile with extended red wing similar to observed ones. Anyway our results indicate that 
magnetic field of CTTSs is signifificantly non-dipole near stellar surface in agreement with 
direct magnetic field measurements (Johns-KruU, 2007). 

4 Conclusion 

We demonstrated that observed intensity and profiles of C IV 1550 doublet lines signifi- 
cantly differ from theoretical predictions based on the assumption that magnetic field of 
CTTSs near stellar surface is close to dipole. We conclude therefore that geometry of 
CTTS's magnetic field near stellar surface is strongly non-dipole. Multipole components 
of global magnetic field of young star or/and small-scale magnetic fields of active regions 
probably produce large divergency of accreted gas stream lines within accretion zone that 
presumbly can explain disagreement between the theory and observations. 

We thank the LOC of the Symposium for the invitation, financial support and hospi- 
tality. 
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